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A circularly polarized light can induce a dissipationless dc current in a quantum nanoring which 
is responsible for a resonant helicity-driven contribution to magnetic moment. This current is not 
suppressed by thermal averaging despite its quantum nature. We refer to this phenomenon as the 
quantum resonant inverse Faraday effect. For weak electromagnetic held, when the characteristic 
coupling energy is small compared to the energy level spacing, we predict narrow resonances in the 
circulating current and, consequently, in the magnetic moment of the ring. For strong fields, the 
resonances merge into a wide peak with a width determined by the spectral curvature. We further 
demonstrate that weak short-range disorder splits the resonances and induces additional particularly 
sharp and high resonant peaks in dc current and magnetization. In contrast, long-range disorder 
leads to a chaotic behavior of the system in the vicinity of the separatrix that divides the phase 
space of the system into regions with dynamically localized and delocalized states. 

PACS numbers: 78.20.Ls, 78.67.-n, 73.23.-b, 75.75.-c 


I. INTRODUCTION 

Nanodevices based on quantum dots, quantum wires 
and quantum rings continue to attract considerable 
attention^— From the physics point of view such sys¬ 
tems are often determined by the interplay of quantum 
interference and charge quantization effects which both 
become more prominent with decreasing system size and 
temperature. Research on electronic phenomena such as 
the Aharonov-Bohm effect, Anderson localization, Kondo 
effect, or Coulomb blockade has been dominating the 
field in the last two decades >^ ~ 15 i 18 ~ — In recent years, 
however, there have appeared numerous proposals to uti¬ 
lize nanodevices in optoelectronics and spintronicsilF— 
This development calls for better understanding of light- 
matter interaction in such systems as quantum wire an¬ 
tennas, artificial atoms, and nanorings. 4 -~— 

One of the main goals of optoelectronics is to design 
and fabricate tunable electronic nanodevices that are ca¬ 
pable of operating in a frequency range unaccessible for 
conventional electronic technologies, i.e. in the so-called 
terahertz (THz) gap. It is widely believed that the fre¬ 
quency gap can be closed using optoelectronic and plas- 
monic devices. There is, however, a serious obstacle for 
such development. The coupling of THz electromagnetic 
field to a single nanosystem appears to be too weak be¬ 
cause the typical dimension of a nanosystem is two or 
more orders of magnitude smaller than the THz wave¬ 
length. A promising way to increase the coupling is to 
use periodic structures (arrays of nanoparticles, grating 
gate structures, multigate structures, etc). Another diffi¬ 
culty originates in dc photoresponse that is only possible 
in the presence of a system asymmetry (which would de¬ 
fine the direction of the photoinduced dc current). In 
two-dimensional systems such an asymmetry might be 
created by boundary conditions^ or induced by a ratchet 
effect (see Ref. for review) . The latter implies a spe¬ 


cial type of grating-gate couplers that could provide the 
required asymmetry. 

Interestingly, the symmetry conditions for photore¬ 
sponse are more relaxed in the multiconnected structures 
such as quantum rings. In particular, the dc circular cur¬ 
rent can be excited in a quantum ring by a circularly 
polarized optical field: £ = £ u exp(— iuit) + h.c., where 
£ is the electric-field component of the electromagnetic 
wave. Such response can be characterized by an orbital 
magnetic moment of the ring 

Moci£ ul x£* u , (1) 

the effect which is commonly referred to as the inverse 
Faraday effect In contrast to other photomagnetic 
effects, the inverse Faraday effect does not involve ab¬ 
sorption of photons or heating, which makes it particu¬ 
larly useful for spintronic applications such as data stor¬ 
age technologies^ Although the magnetic moment gen¬ 
erated in a single ring is relatively small, an ensemble of 
nearly identical quantum rings may give rise to large opti¬ 
cally controlled macroscopic magnetization. Two points 
are especially important in view of possible applications: 
(i) the proportionality coefficient in Eq. © is an odd 
function of frequency, so that the effect is sensitive to 
the helicity of polarization, and (ii) the effect is sizable 
even in the limit of long wavelength such that £ bJ does 
not vary within the ring dimension. Hence, quantum 
nanorings and ring-based arrays can be used as effective 
helicity-driven sensors for THz radiation. 

Historically, the inverse Faraday effect has been pre¬ 
dicted by PitaevskiiSi and first observed by van der Ziel 
et al£2. Much of the current interest to the phenomenon 
originates, however, in the experiments by Kimel et 
al£2r— on ultrafast femtosecond magnetization dynam¬ 
ics in thin ferrimagnets. In this paper we leave aside 
many unresolved issues in the theory of the inverse Fara¬ 
day effect in magnetic materials but focus instead on the 
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excitation of current and magnetic moment in a single¬ 
channel quantum ring. This problem has been analyzed 
recently by Kibis^ using perturbative analysis (see also a 
more recent publication^) and by Alexeev et using 
a master equation while disregarding diamagnetic cur¬ 
rent. A similar system but in the presence of strong 
spin-orbit interaction at zero temperature has been re¬ 
cently considered— The electric dipole moment oscil¬ 
lations in quantum rings at a finite temperature^ and 
inverse Faraday effect due to the flux change through 
classical (macroscopic) metallic rings subjected to short 
optical pulses were also discussed i 31 ’ 32 The inverse Fara¬ 
day effect in mesoscopic chaotic cavities has been stud¬ 
ied in Ref. Is^L In this paper we focus on the inverse 
Faraday effect in quantum rings. In contrast to previ¬ 
ous publications we develop a nonperturbative approach 
that remains valid for the case of strong coupling to elec¬ 
tromagnetic field. We focus specifically on the resonant 
enhancement of the inverse Faraday effect in the absence 
of spin-orbit interaction and at relatively high tempera¬ 
tures. 

Optically-induced circular current / ra( j has a number 
of similarities to persistent current / per . The latter may 
flow in a quantum ring at thermodynamic equilibrium. 
Both currents are dissipationless and vary periodically 
with magnetic flux piercing the ring. Both currents 
arise due to time-reversal symmetry breaking by mag¬ 
netic field and/or by circularly polarized light. Conse¬ 
quently, / per is an odd function of the magnetic field 
Ipei{<P) = while / ra d changes sign upon the 

inversion of both magnetic field and optical field helicity 

Aad(/*>kO = ^rad( ^)- (2) 

Here we introduce </ = ‘F/'Fo, where $ is the magnetic 
flux piercing the ring and $0 = he/e is the flux quantum. 

The persistent and optically induced currents are, how¬ 
ever, completely different when it concerns their temper¬ 
ature dependence. The averaged persistent current is 
exponentially suppressed with increasing temperature T 
above the level spacing at the Fermi level Anamely 
/per oc exp(— T/Ap), while the optically induced cur¬ 
rent decays much slower, / ra d oc Ap/T, as we show 
below. Another closely related difference is related to 
the role of mesoscopic fluctuations in these two currents. 
Such fluctuations provide a dominant contribution to the 
persistent current for all temperatures, which makes it 
sensitive to the type of the thermodynamic statistical 
ensemble— (Since fluctuations exponentially exceed the 
averaged value of the current, the quantity / per is not 
representative for a given isolated ring.) 

In contrast, as we demonstrate below, the mesoscopic 
fluctuations of / ra d are small for temperatures exceeding 
the mean level spacing so that the ensemble-averaged op¬ 
tical current is well defined. The dependence of / ra d on 
the type of thermodynamic averaging is, therefore, neg¬ 
ligible at high temperatures. 

To conclude the comparison of / per and / ra d, we note 
that persistent current might show up indirectly at high 


temperatures. In particular, it was demonstrated in 
a series of publications^ - — that the tunneling current 
through a single-channel quantum ring is blocked by per¬ 
sistent current. This effect, caused by an interplay of 
quantum interference and charge quantization, has been 
named the persistent-current blockade (PCB) in analogy 
with the well-known Coulomb blockade. In contrast to 
the latter, the PCB persists for much higher tempera¬ 
tures despite its essentially quantum nature. The meso¬ 
scopic fluctuations of / per in the regime of PCB lead to 
the splitting of Aharonov-Bohm resonances at high tem¬ 
peratures. Similarly, we will find that the current / ra d 
survives up to sufficiently large temperatures. 

The slow decay of the inverse Faraday effect with tem¬ 
perature yields additional advantages for optoelectronics 
and spintronics. Therefore, in this paper we focus on the 
high-temperature regime 

T> A f . (3) 

We calculate / ra d for arbitrary coupling to electromag¬ 
netic radiation by paying a particular attention to res¬ 
onance effects. In the weak field limit we predict se¬ 
ries of narrow resonances in the frequency dependence of 
/ ra d- Each resonance corresponds to excitation frequency 
coinciding with the distance between neighboring levels. 
Weak short-range disorder splits the resonances and in¬ 
duces particularly sharp resonant peaks in magnetiza¬ 
tion. For the case of large field we find using the quasi- 
classical approximation that the resonances broaden and 
merge into a single wide peak. The width of the peak is 
limited by a thermal band for moderate coupling while 
it is proportional to the square root of the wave ampli¬ 
tude for very strong fields. In a clean limit, i.e., in the 
absence of disorder, the corresponding circular current is 
dissipationless. The presence of long-range disorder leads 
to a chaotic behavior of the system in the vicinity of the 
separatrix that divides the phase space into regions with 
dynamically localized and delocalized states. In contrast, 
weak short range disorder leads to the appearance of ad¬ 
ditional particularly sharp and high resonant peaks in dc 
current and magnetization. 


II. IDEAL QUANTUM RING 

An ideal quantum ring placed in a circularly polarized 
electromagnetic field can be described with an effective 
stationary Schrodinger equation by transforming to the 
rotating frame. The solution to this equation is straight¬ 
forwardly obtained in two limiting cases: (i) for weak 
electromagnetic field, such that the coupling energy is 
small compared to the level spacing, and (ii) for strong 
field, such that the level spacing in the ring is negligible 
compared to the coupling matrix element between elec¬ 
trons and photons. 

Before going into details let us recall first the well- 
known concept of the persistent current in a ballistic 


3 



u 


FIG. 1: Single-channel quantum ring threaded by magnetic 
flux Circularly polarized light induces a nondissipative 
current circulating around the ring. 


single-channel ring. In the absence of both electromag¬ 
netic radiation and magnetic flux each energy level in an 
ideal ring carries the electric current I n = Iqu, where 


Io = 


eh 

2ttMR 2 ' 


(4) 


acting on the electrons in the ring is homogeneous and is 
given by 


£u « £o(e x - ie y )/2, (7) 


where e x and e y are unit vectors in x and y directions, 
respectively, and £q is the amplitude of the field. The 
Schrodinger equation for the ring, which is threaded by 
a magnetic flux, is given by 

n,Tf /pi \ 2 

= -y (- i<j>) f - e£ 0 R cos (p - wf)'F, (8) 


where £o = hr /MR 2 and p is the polar angle shown in 
Fig. [0 In Eq. © we neglect small corrections arising 
due to a finite size of electron wave function in the radial 
direction. The function corresponds to a state which 
carriers a dc current given by 


I = 2t r/ 0 


2 i \ dp dp J 



(9) 


Here R stands for the radius of the ring, M is an effective 
electron mass, and integer number n numerates energy 
levels. Due to the evident symmetry I n = /_„ the total 
equilibrium current circulating in the ring vanishes. If 
the ring is threaded by a magnetic flux one finds 

I n = Io(n-4>)- (5) 


where (•••}( stands for time averaging. 

Equation (|5]) can be transformed into a stationary 
Schrodinger equation using a rotating reference frame, 

(<p, t) = e -iEt e i<t, Ve in^- U t) x ^ _ ( 10 ) 

where we introduce the dimensionless frequency and cou¬ 
pling 


Thus, for generic flux the exact cancellation is absent 
and a dissipationless persistent current flows. At zero 
temperature one estimates the persistent current as 

Iper = 'y ( In ^ 0 , ( 6 ) 

where En' 1 are energy levels in the absence of radiation 
and Ep is the Fermi energy. It is evident from Eq. (0) 
that the persistent current flows even in the absence of 
an external electric or electromagnetic field. 

The persistent current (0) has been indeed observed in 
experiments with ensembles of nanorings.— (see also 
Ref. Fill for review). Quantitative theoretical explanation 
of these experiments is, however, much more involved 
given that the rings are typically disordered and not one¬ 
dimensional while the effects of electron-electron interac¬ 
tions are not negligible. 

Nevertheless, we shall start with the discussion of the 
simplest model, which is a clean single-channel quantum 
ring, and postpone the generalization of our results for 
the disordered case to Sec. m The case of a multichannel 
ring as well as the effect of electron-electron interaction 
will be discussed qualitatively in Secs. IVIland IVIII 

We consider a single-channel nanoring subject to a cir¬ 
cularly polarized radiation with frequency to. The radius 
of the ring is naturally assumed to be small compared 
to the wavelength of light. In this case the electric field 


oj e£o R 

Tl(jj — 5 OL — 

£o £o 


(ii) 


Here and in what follows we put h = 1. The eigenenergy 
is conveniently parameterized by 


E = 


£ T cr — 



£ 0 - 


( 12 ) 


where e is a dimensionless energy. 

The wave function in the rotating reference frame, 
x(0), obeys the differential equation 

X" + 2x(£ -W) = 0, (13) 


where the double prime stands for the second derivative 
with respect to the angle 9 = <p — uit , while the effective 
potential is given by 

W{6) = — a(l + cos 9). (14) 

This potential is plotted in Fig. [2] The boundary condi¬ 
tions for Eq. 03 read 

x(0)/x(2tt) = x'(0)/x'(2 7 r) = e 2 ^+""). (15) 


Thus, in the case of an ideal quantum ring the problem 
is reduced to the solution of the Schrodinger equation 
which corresponds to a quantum physical pendulum with 
nonperiodic boundary conditions. 
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the laboratory frame, 


$(0) = e -iE^t e 


imp 




£ ,„ ) = £ (18) 


The current corresponding to n —th level in the absence 
of radiation is given by Eq. ([5]). By weighting these cur¬ 


rents with the corresponding Fermi function fp 



we 


arrive with the help of the Poisson summation formula 
at the expression for persistent current in a clean ring, 


FIG. 2: Potential of quantum pendulum. 



FIG. 3: Energy levels in laboratory and rotation frames. 


The solution to Eq. m supplemented with bound¬ 
ary conditions of Eq. (1151) gives rise to the eigenenergies 
e n and eigenfunctions Xn- The corresponding radiation- 
dressed functions 4' n (</j, t) are, then, found from Eq. (fTOli . 
where E = E n is related to e n by Eq. m These func¬ 
tions represent a full basis for electron states in the quan¬ 
tum ring. The energy levels in both laboratory and ro¬ 
tated frame are illustrated in Fig. [3] 

The wave-function \I> n is also an eigenstate of the cur¬ 
rent operator which corresponds to the dc current ex¬ 
pressed in terms of \ n as 


I n = h dd 


■ d\n 

1 86 


9Xn \ , |2 

Xn qq ] + n u\Xi 


. (16) 

In the absence of radiation, i.e.,for a = 0, we simply 
obtain 


Xn 


E n 


gi(n—ncj—cf>)0 

e 0 {n-(j)) 2 

2 


£n 



win — (f>). 




(17a) 

(17b) 


From Eqs. m it is easy to see that the conventional 
results for a ballistic single-channel ring are restored in 


I P ei(<j>) = ^2 ImSm(2Trm(j)), (19) 

771— — OO 

where I m = —Jo f dx x sin(2nmx) fp(£oX 2 /2). One can 
see that J pe r(<)>) decreases exponentially with tempera¬ 
ture [proportional to exp(— T/Ap )] and J per (0) = 0. The 
total current I = J per + J ra d, however, includes an addi¬ 
tional radiation-induced contribution J ra( j. 


III. RADIATION-INDUCED CURRENT IN THE 
WEAK COUPLING REGIME 

Let us now turn to the radiation-induced contribution 
to the current in the clean ring in the case of relatively 
large temperatures such that 

£o A p CT< Ep , (20) 

where A^ ~ Eonp is the level spacing at the Fermi level 
and tif — (2Ep /eq) 1 / 2 ^$> 1. In contrast to the persistent 
current, the radiation-induced contribution J ra( j is not 
exponentially suppressed in this regime. Still, similarly 
to persistent current, J ra d varies periodically with <j> and, 
therefore, can be tuned by external field. 

We consider first the case of an isolated ring disre¬ 
garding coupling to the thermal bath. The amplitude of 
radiation is also assumed to be switched on adiabatically. 
Later on, we generalize the obtained result to account for 
relaxation processes. 


A. Isolated ring, adiabatic radiation switching 

For adiabatic switching a one-to-one correspondence 
between unperturbed quantum states and radiation- 
dressed eigenfunctions can be established. Namely, the 
states described by in Eq. m transform adiabati¬ 
cally into Assuming naturally that the unperturbed 
system was in a thermal equilibrium we arrive at the fol¬ 
lowing result: 


OO 

Jrad = ]T SI n (21) 

n=—oo 


where 5I n = I n — In , I n is the current corresponding 
to radiation dressed functions 4' n , and f n = fp[En ^ is 
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the equilibrium distribution function over unperturbed 
energy levels. The result of Eq. m can be rewritten as 

OO 

/rad= Jn(fn~fn+ 1) , (22) 

n=— oo 

where we introduced 

n 

Jn= Y, SIm ' (23) 

m =—oo 

In the limit of relatively large temperatures, such that 
inequalities ([20]) hold, one can further simplify Eq. ([22]) 
as 


/rad — ^ ^ Jr 


n=— oo 
oo 


dfn 

dn 


(24) 


- E 


Jn 


4n* cosh [(n — n F — (j))/2n*] 


where 


n* = T/A f > 1 (25) 


is the number of quantum levels in the temperature win¬ 
dow. The term <f>/2n* in the argument of hyperbolic 
cosine is small but it is needed to preserve the exact in¬ 
variance of Eq. (l24ll with respect to a shift of magnetic 
flux by a flux quantum: (f> —> <j) + 1. 

Let us now evaluate J n for the case of weak coupling to 
external radiation, i.e., for a -C 1, by taking advantage 
of perturbation theory. Keeping terms up to the second 
order with respect to a one finds the spectrum 


{n-n u - 4>) 2 
S n — --- & + 


a 


4 [(n -n u - cj)) 2 - 1/4] 


(26) 


in accordance with earlier work by Kibis^ We note that 
optical field induces a change in the current for each 
radiation-dressed quantum level (this effect was not dis¬ 
cussed in Ref. |26j). To the second order in a [with the 
same precision as in Eq. (l26|) ] we obtain 


Sin = ~h 


a 


n — n u - < 


2 [(n -n u - (j)) 2 - 1/4] 2 
Substitution of Eq. (071) into Eq. (1551) yields 


(27) 


Jn — Tn 


4 (n — n u - (j) + 1/2) 


a 

2 - 7x2 > 


= Io 


772. (28) 


where we introduced 

S n = — -— = Hu — n + <f> — 1/2, (29a) 

£o 

A n = - E& = so (n - + 1/2). (29b) 

The energy A„ is nothing but the spacing between the 
level n + 1 and n. 


It is evident from Eq. (l28l) that J n is strongly enhanced 
provided a resonance condition u> ss A„ for a given n. Let 
us assume that such a resonance takes place for n = N 
such that Sn -C 1. In this case one also finds |5„| > 1 
for all n ^ N. Consequently, the current is dominated 
by the contribution coming from transitions between the 
levels N and N + 1. It is immediately concluded from 
Eq. 051) that the perturbation theory applies for a <C <5, 
but fails in the opposite limit. Let us, therefore, modify 
Eq. 051) to take into account non-perturbative effects. 

In order to evaluate the resonant contribution to the 
current let us for a moment neglect all optically-induced 
transitions except for the transition between levels N and 
N+1. In the rotating wave approximation, the corre¬ 
sponding two-level Hamiltonian reads 


£jV W 
W £n+1 


(30) 


where the wave functions in the rotation frame \n, 'Xn+i 
and the corresponding energies £n+n+i are given by 
Eqs. (I17al) . while W stands for the matrix element of the 
optical transition TV + 1 o IV, 


P2TT 

W = — a / e lB cos (6)dd = — —. (31) 

Jo 2 

The eigenfunctions of the projected Hamiltonian (1301) are 
given by 


Xn = 


Xn+i = 


XN - 13 XN +1 


VT+P' 


Xn+i 


•2 

P Xn 


VlT P 2 


where we introduced yet another parameter 
a sign(<5iv) 


P = 


I^jvI + V Sn a 


(32a) 

(32b) 


(33) 


The phase factors in Eqs. (l32l) are taken in such a way 
that functions \n and xat+i transform, respectively, into 
Xn and Xn+i for both positive and negative Sn for a —> 
0 . This leads to appearance of a modulus of |< 5 jv | and 
sign(dAf) in Eq. (l33l) . From Eq. (flfil) one obtains the 
result for currents 


SIn = —SIn +i = Jo t '' o 2 ’ ( 34 ) 

1 + P z 

which is illustrated schematically in Fig. [I] One can see 
from Eq. (1341) that the current variations SIn and SIn +i 
can be as large as Iq. On the other hand the contribu¬ 
tion of other transitions with n ^ TV is suppressed by a 
small factor a 2 and can be estimated as ck 2 To- Thus, the 
dominant contribution to Eq. (1221) indeed comes from the 
transition with n = TV. From Eq. (1341) we find 


Jn — SIn — 


To_ ^ _ 

2 yja 2 + <5^(|5jv| + V^+S%) 


(35) 
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Following the procedure described in the beginning of the 
section we neglect transitions between levels that vary 
adiabatically and, therefore, refer to the equilibrium dis¬ 
tribution of the unperturbed system (i.e.,to the state be¬ 
fore the external radiation is adiabatically switched on). 
By doing so we arrive at the following expression for the 
radiation-induced current, J ra( j ~ Jjv (/at — /jv+i), for 
to ~ A 


equality in the thermal factor in Eq. (1371) one finds 


Tenv (^) 


8 n* cosh 2 [(n w — uf — l/2)/2n*] 


rmax 

-'rad 


cosh 2 [(oc — Af)/2Sw] 


(39) 


Trad — -o-, (36) 

4n* cosh 2 [(N — uf — <f>)/2n*] 

where we take into account that in the adiabatical case 
the current SIn must be weighted with the unperturbed 
Fermi distribution function /„. 

Summing up the contributions from all levels we ar¬ 
rive at a more general result which includes both non¬ 
resonant and resonant contributions 


Trad — 


X 


lO ^ a 

w+7^+1 

i 

4n* cosh 2 [[n — uf — (f>)/2n*] 


(37) 


The dependence of the current / ra d given by Eq. (1371) on 
frequency is shown in the upper panel of Fig. [5] Posi¬ 
tions of the peaks are found from the conditions S n = 0. 
The smooth envelope of the peaks is due to the ther¬ 
mal factor l/4n* cosh 2 [(n — uf — <j>)/2n*]. The central 
peak corresponds to resonance excitation of levels at the 
Fermi energy while its amplitude corresponds to a max¬ 
imal possible optical response for weak coupling, which 
can be estimated as 


[max 

-'rad 


-fr) 

8 n* 


IqAf 
8 T 


(38) 


The distance between resonant peaks is given by £o- Re¬ 
markably, the coupling strength a drops out from the 
result of Eq. (1351) . Thus, the radiation-induced con¬ 
tribution to current might be large even in the weak- 
coupling regime, and decays as T -1 in contrast to the 
persistent current contribution which decays exponen¬ 
tially with temperature. 

It is worth stressing that positions of resonances de¬ 
pend on <j> and are therefore tunable by magnetic flux 
piercing the ring. Changing the magnetic flux by the 
flux quantum, (f> —> <f> + 1 is equivalent to the substitu¬ 
tion S n —> 8 n -h which proves that the result of Eq. © 
is a periodic function of </> with the period 1. Thus, in¬ 
stead of varying frequency of radiation one can probe 
optically-driven resonances in the ring by varying exter¬ 
nal magnetic field. 

Let us present analytical expression for the smooth en¬ 
velope /env(w) of the resonance peaks plotted with the 
dashed line in the upper panel in Fig. [5] The maximal 
peak values are found from the condition S n = 0, which 
is equivalent to n « n u + (f> — 1/2. Substituting the latter 


where 8to = eon* = T/uf is the envelope width. From 
the physics point of view the frequency range Sui corre¬ 
sponds to the variation of the level spacing within the 
temperature window, hence Sto — (dA n /dn)n*. 

It is worth noting that Fig. [5] corresponds to the case 
of positive hclicity (to > 0). The current changes sign un¬ 
der simultaneous inversion of helicity and magnetic field 
[see Eq. ©], he. ,under the time reversion. Consequently, 
the field-independent contribution to the current changes 
sign under the inversion of helicity only. 

Before closing this subsection, we shall briefly discuss 
the role of mesoscopic fluctuations. As we mentioned 
above, such fluctuations dominate persistent-current es¬ 
pecially at high temperatures, when the averaged value 
of /p er is exponentially small. Let us estimate the fluc¬ 
tuations of J ra d using simple arguments. First of all we 
shall notice that only the electron levels within the tem¬ 
perature window around the Fermi level can be popu¬ 
lated or depopulated. The number of such levels is es¬ 
timated as n* [see Eq. Thus, the fluctuations of 

the total number of electrons in the ring are given by 
AN ~ yfnF ~ \/T/Ap- The corresponding fluctuation 
of the chemical potential reads A fj, = ApSN ~ y/TA f■ 
Such fluctuations do not affect our results for / per pro¬ 
vided the factor df n /dn entering Eq. (l24l) does not fluc¬ 
tuate much. This is indeed the case for A fi/T <C 1. In 
this limit we estimate the mesoscopic fluctuation of the 
optically-induced current as 


A / ra d A/li / A F 1 

-/rad - — ~T~ ~ V ~T~ « 1 


(40) 


Thus, we conclude that for high temperatures meso¬ 
scopic fluctuations of / ra( j are suppressed.— This, in turn, 
implies that, in contrast to persistent current, / ra( j is not 
very sensitive to the choice of thermodynamic statisti¬ 
cal ensemble. The robustness of / rac j with respect to the 
mesoscopic fluctuations can be understood rather easily 
from a simple physics argument. Indeed, the resonant 
condition for a pair of neighboring levels in the tempera¬ 
ture window near the Fermi surface is satisfied provided 
that interlevcl distance equals the radiation frequency. 
Since the distance between levels does not depend on 
ensemble, the only condition required is that the fluctu¬ 
ations do not move the majority of relevant levels out of 
the temperature window. This condition is equivalent to 
T > A/i. 
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FIG. 4: Currents of quantum levels for resonant excitation 
(w « An) and weak coupling to radiation (a< 1). 



FIG. 5: Upper panel: Radiation-induced circulating current 
for the case of weak coupling to radiation. The resonances 
corresponding to excitations of different pairs of levels are 
well separated; Lower panel: Resonances overlap and merge 
into a single wide peak for the case of strong coupling. 



FIG. 6: Currents of quantum levels for strong coupling to 
radiation (a 1). 


system as 

Pn = -ct£ 0 lm[e %ut p] + 7 (f N - p N ), (41a) 

Pn+i = a£ 0 Im.[e lut p) + 'y(f N+1 -p N+1 ), (41b) 

-P = (*A N + 7 v )p + ^(Piv +1 - PN)e~^\ (41c) 

where pn+i = Pn+i,n+i, Pn = Pn,n, and p = Pn+i,n 
parametrize the relevant components of the density ma¬ 
trix, while 7 and 7 ^ are the relaxation and dephasing 
rates, respectively. In what follows we substitute 

p = e~ iut p (42) 


and search for stationary solutions of Eqs. m■ Straight¬ 
forward analysis yields 


SpN 


P 


a 2 r y (/jv + i — /at) 

2[T(S%+Tl) + a 2 T lp } 
ol(pn+i — Pn) 

2(S N + iT v ) ’ 


-SpN+i, (43a) 
(43b) 


B. Ring coupled to thermal bath 


In this subsection we shall turn to the case of ideal 
quantum ring coupled in addition to a thermal bath while 
still assuming a -C 1 , i.e. a weak coupling to electromag¬ 
netic radiation. We will see that in this case qualitative 
physical picture drown in Fig. [5] and in previous subsec¬ 
tion remains intact provided that the rate of relaxation 
is sufficiently low. The only essential difference is related 
to the shape of the optically-induced resonances in the 
current. 

Similarly to the previous section we assume <5 tv ■C 1 
and consider the only pair of resonant levels N and N+l. 
In this regime the effect of the thermal bath can be de¬ 
scribed by the simplest Markovian model, which is ex¬ 
pressed in terms of the density matrix for the two-level 


where Sp m = Pm — f n is the radiation-induced variation 
of level population while T = y/eg and T^, = y^/eg stand 
for dimensionless rates. The corresponding expression for 
dc current reads 

OO 

4 r ad =Io ^2 5 Pn,n- (44) 

n =—oo 


Substituting Eq C3l) into Eq. (El) , we obtain the result 


Aad — 


/q a 2 ry(/Af - f N + l) 

2 r(6% + T£) + 


(45) 


Thus, the resonant radiation-induced current in the pres¬ 
ence of a thermal bath is still given by Eq. (l36l) with 


Jjv = 


crY a 


2 T(5 2 N +Tl) + a 2 T v 


(46) 






















The inequality > T/2 holds since relaxation of level 
populations leads to dephasing. Assuming that there are 
no other sources of dephasing one can simplify Eq. (14611 
to 

Jjv = T^ + rV4 + aV2- (47) 

By comparing Eq. (H71) with Eq. (1551) we conclude that 
the maximal value of Jn is smaller for the ring cou¬ 
pled to a thermal bath than for an isolated ring by 
a factor a 2 /(a 2 + T 2 /2) as intuitively expected. That 
explains why the resonance width increases from a to 

yj a 2 + T 2 / 2 . 

For small relaxation rate, T<Cl, the dc current shows 
a series of sharp resonances in full analogy to the case 
of isolated adiabatic system. In the limit T —> 0 the 
only difference between these two cases is related to the 
different shapes of resonances. For a finite T the differ¬ 
ence becomes more essential due to the presence of some 
dissipation caused by interaction with the bath. The dis¬ 
sipated power P can be estimated from the conventional 
formula for Joule heating 

P = R J ( 48 ) 

where £ v (y, t) = —£q sin(c/j — cut) is the projection of the 
electric field on the current direction, R is the ring radius, 
and the current is given by a generalization of Eq. j9]) for 
the system described with the help of the density matrix, 

/(¥>,*) = IoJ2pnm(t) - A ( 49 ) 

n,m ' ' 

In contrast to Eq. m we, however, should keep in 
Eq. (l49l) the ac contribution to the current. It is pre¬ 
cisely this contribution that is responsible for dissipation. 
Straightforward analysis yields the dissipated power 


We recall that resonances obtained in the weak- 
coupling regime (see upper panel of Fig. 0 are separated 
by the distance £o and have the width ae o. Therefore, 
the resonances are expected to overlap if coupling to radi¬ 
ation increases. In the regime of overlapping resonances 
the two-level approximation used in the previous section 
is no longer justified and a more accurate analysis has 
to be performed. Simple consideration below shows, in¬ 
deed, that for a 1 the dependence of radiation-induced 
current on frequency is given by a single peak of a large 
amplitude. This dependence is depicted schematically in 
the lower panel of Fig. 0 

For an isolated ring strongly coupled to radiation the 
effective potential W{6) is large enough to localize the 
states near 9 = 0 in the rotating frame. These localized 
states correspond to the energy range e < 0 in Fig. [2] 
Simple quasiclassical analysis of Eq. (fl3l) shows that the 
total number of localized states is of the order of y/a while 
the distance between the levels is proportional to y/aeo- 
These energies are close to the bottom of the parabola 
in Fig. [2] and correspond to n lying in a vicinity of n u 
such that \n — n u — <j)\ < y/a. In the laboratory frame the 
localized states form a band of the width y/a centered 
around n = n^. In the absence of disorder all states 
in the band have a certain chirality. The helicity of the 
radiation determines the sign of n u and, consequently, 
the chirality of the localized band. 

Since wave functions Xn(9) for localized states are real 
the only term which contributes to the radiation-induced 
current is the last one in Eq. ( 110 . Using the normaliza¬ 
tion condition J \x n {0)\ 2 d6 = 1 we estimate the localized 
state contributions to the current I n ps = constant 
for \n — riuj — <j>\ < y/a\ as illustrated in Fig. 0 On the 
other hand for energies outside the localized band, i.e.,for 
\n — riu; — <j>\ y/a , the radiation does not affect the cur¬ 
rent in any essential way hence the perturbative result 
of Eq. (1271) remains valid. With the help of Eq. (E3l) we 
obtain 


P = 2ttR£qIo 


aT (N + l/2-<f>)(f N -f N+1 ) 
4 S 2 N +T 2 /A + a 2 /2 


(50) 


which is proportional to both the coupling constant a 
and the rate F which characterizes the coupling to the 
bath. 


IV. STRONG COUPLING TO THE RADIATION 
A. Isolated ring, adiabatic radiation switching 

In this section we focus on the strong-coupling regime 
a ^$> 1. Similarly to the previous section we consider 
first the case of a completely isolated ring assuming that 
radiation switches on adiabatically so that we can deduce 
the level occupation numbers from an equilibrium state 
at an initial moment of time. We, then, turn to the case 
of a ring coupled to a thermal bath for which we do not 
need to make such an assumption. 


Jn 


4 (n—n^ —(f>)‘ 2 5 

Cia - , 


|n-n w - 0 | > y/a, 

(51) 

|n—riuj — 0| <C y/a, 


where Ci ~ 1 is a numerical coefficient. Let us now as¬ 
sume that the temperature is sufficiently large so that 
n* y/a, or, equivalently, T » y/aAp. Then, the dis¬ 
tribution function does not change within the width of 
the band. Substituting Eq. (l5lT) into Eq. (l24l) and replac¬ 
ing the summation over n with integration we estimate 
the current as 


trad 


/ 0 a 3/2 

n* cosh 2 [(noj — n_p)/2n*] 


(52) 


This result suggests that the maximal value of the current 


max 

-'rad 


W /2 

n* 


(53) 


is achieved for u> ps Ap while the width of the broadened 
resonant peak is proportional to T/tif. Comparing this 
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result with that of Eq. (1551) we conclude that the maximal 
achievable current increases with the radiation strength. 

For sufficiently large values of a , such that T -C 
the dependence of the derivative df n /dn on n 
becomes stronger than that of J n . In this case the cur¬ 
rent is given by / ra d ~ J np , where J n is still determined 
from Eq. The maximal value of the current in this 

case reads 


/™d X ~ 7 0 a, for riF»v / a>n*- (54) 

In the limit of very large coupling y/a n Fl the maximal 

current saturates at the value 

-Cd X ~ hn 2 F , for y/a^>n F , (55) 


where higher harmonics with respect to the angle 6 are 
neglected. Substituting Eq. (1551) into Eq. (1551) we find 




J n +1 


OL£o Fu-\-l F n 

2 S n 


(59) 


where the function F n is determined by the balance equa¬ 
tion 


a 2 T v / F n+ i — F n F n _i~F n \ 
2 ^ <5 2 +r 2 + 5l_ 1+ T%) 


r(F n -/„), ( 60 ) 


with the same definition of dimensionless collision rates 
T = y/eo and = 7 ^/£o as in the previous section. 
Since both collision rates are large T 1, T v 1 
Eq. (HOD can be rewritten in the differential form as 


while the width of the peak in both regimes (1551) and 
(1551) is given by ydaeo. Equation (1551) has a very clear 
physical sense. For such a large coupling all electrons 
are localized in the rotation frame and rotate with the 
velocity v F cx n F . The current is given by the ratio of 
the total charge Q ~ en F to the time of the electron 
traveling around the ring, which is proportional to l/v F . 
This yields Eq. (1551) . 


B. Ring coupled to thermal bath 


8 ( 1 dF n \ - t 

2 T dn V<5n +T'f dn ) n fn 


(61) 


The solution to Eq. ED can be found using the ansatz 


F n = fn + T^Gn^-, (62) 

which is justified for sufficiently high temperatures such 
that n* T v . Substituting Eq. (1551) into Eq. (iGll) and 
neglecting terms which are proportional to d 2 f n /dn 2 and 
d 3 fn/dn 3 we arrive at the following equation for G = G n : 


Similarly to the previous section the radiation-induced 
current in the presence of the coupling to a thermal bath 
can be calculated using Eq. (1551) . The two-level approx¬ 
imation used above is, however, no longer justified due 
to a strong overlap of resonances corresponding to opti¬ 
cal transitions between particular energy levels. In this 
regime the relaxation of diagonal and off-diagonal ele¬ 
ments of the density matrix is governed respectively by 
the terms 7 (/„ — p n ,n) and — 7 v p n , m in the collision inte¬ 
gral. For simplicity we shall neglect the possible depen¬ 
dence of collision rates 7 and r y v on energy. 

The equation on the density matrix takes the form 


dF n 

dt 


= [(n 


(j))e 0 - w] 


dF n 

d6 


ie 0 d 2 F n 
2 dd 2 


+ (. F n+l e w + F n _ ie ~ ie - 2 F n cos9) 


+ 7 (fn - F n ) + 7 v{F n - F n ), (56) 


where F n = p n , n = 


± ST ddF n{0) and 


F n (t, e) = J2 Pn,n+ke ik ( e - ut) . (57) 

k 

The Equation (1551) is easily analyzed provided 7 > eoj 
7 V eo- In this limit we search for a stationary solution 
to Eq. (1551) in the following form 


Fn — F n + a. n e l6 + (3 n e 16 , 


VG 


d \ 1 / dGY 

dx 1 + x 2 \ dx J 


(63) 


where x = 5 n /T v and 77 = 2rr^/a 2 . The dimensionless 
parameter 77 characterizes the strength of thermalization 
rates relative to optical transition rate and can, therefore, 
be regarded as a measure of thermalization intensity. For 
77 1 one can neglect the term dG/dx in the right hand 

side of Eq. (1551) . Thus, for relatively fast thermalization 
we obtain 

2 t 

Ga %(i + =er tor,,>>1 - <64) 

This result also applies for |x| !2> rj ~irrespective of 
the value of 77 . For \x\ -C 77” 1 / 4 one can disregard the 
left hand side of Eq. (1551) . Thus, the behavior of G in the 
limit of relatively slow thermalization is given by 


G 



|x| < 77 - 1 / 4 , 
\x\ » 77 - 1 / 4 , 


for 77 -C 1 . 


(65) 


The results of Eqs. (1551) and (1551) have to be substituted 
into Eq. (1621) in order to obtain the distribution function. 
Using the latter in Eq. (ED, taking into account that 
F n = p n ,m and replacing the summation over n with 
integration we arrive at the following result for current 


Trad 


_ Iq_ 

cosh 2 [ 


Tluj —TlF 1 
2 n* \ 


a 2 /r, 

a 3 G(T v /T) 3 /\ 


a 2 « r3 r, 
a 2 » rl r. 
( 66 ) 


(58) 


n 
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It is worth noting that in the limit ~ T and a>T 2 
the second line in Eq. coincides with the result of 
Eq. (l52ll obtained for the adiabatic case. We should also 
note that in the derivation of the asymptotic behavior 
for large a expressed by the second line in Eq. (l66l) the 
large-temperature limit, T Ahas been implicitly 
assumed. 


V. DISORDERED RING 


Surface roughness, impurities and external Coulomb 
potentials are the main sources of disorder which are al¬ 
most impossible to avoid in a realistic nanoring. The 
disorder leads to backscattering that would invalidate the 
analysis of the previous sections. In this paper we do not 
consider the limit of strong disorder such that Anderson 
localisation on the scale of the nanoring circumference 
sets in. Instead we focus on the cleanest possible but 
still realistic systems for which disorder can be regarded 
as small, i.e.,the corresponding mean-free path is large or 
comparable with the ring radius. We further distinguish 
the cases of short-range and long-range disorder. 

The short-range disorder leads to a scattering between 
right- and left- moving electrons. One may naively ex¬ 
pect that such processes would merely lead to addi¬ 
tional broadening of the resonant peaks for the radiation- 
induced current. Contrary to expectations the weak 
short-range disorder is shown below to split the reso¬ 
nances and induce new narrow resonance peaks with the 
amplitude enhanced by a large factor of the order of uf 
as compared to that in the case of an ideal ring. The 
suppression of these new resonant features happens only 
for sufficiently strong disorder. 

The effect of long-range smooth disorder is entirely dif¬ 
ferent. A long-range potential does not affect the result 
for radiation-induced current in the limit of weak cou¬ 
pling a -C 1 , but may affect optical transitions if the 
light amplitude is sufficiently large. In the latter case the 
physics of the system is equivalent to those of a physical 
pendulum described within the quasiclassical approxima¬ 
tion. Random long-range potential leads to a classical 
chaotic behavior of the system, which results in the ap¬ 
pearance of a thin chaotic layer near the separatrix of the 
physical pendulum. 


A. Weak short-range disorder 

Let us start with a more detailed analysis of the model 
in the presence of weak short-range disorder. In an ideal 
ring all energy levels are double degenerate provided the 
dimensionless magnetic flux (j> is an integer or half-integer 
number. In both cases every level with a positive chiral¬ 
ity has a partner with a negative chirality which corre¬ 
sponds to the same energy. A short-range disorder in¬ 
duces backscattering which prevents the use of chirality 


as a quantum number and mixes the pairs of degenerate 
states. 

For the sake of definiteness we shall focus on the vicin¬ 
ity of cj) = 0 first. We further assume that uj « A at and 
aCl. One may still remember from the analysis of the 
previous section that under such conditions the resonant 
optical transition between N —th and (TV + 1)—th level is 
the only one which is relevant in an ideal ring. The pres¬ 
ence of disorder potential, U(<p), mixes the states of pos¬ 
itive and negative chiralities. If both disorder and radi¬ 
ation are sufficiently weak (the corresponding conditions 
will be formulated below) the resonances in radiation- 
induced current can be obtained within a 4-level approx¬ 
imation based on the Hamiltonian projected on the states 
±N and ±(TV + 1), 


H = 


where 


/ 

171 ( 0 ) 

-C/jV 

ae oe lult /2 

U* N 

ae 0 e ~ iuit /2 

rj N +1 

0 


Un 

0 

rj -N 

V 

0 

Un+i 

0 


0 \ 


77* 

U N+1 

0 

e(°) 

-W-fib 


(67) 


Un = 


2tt 


U{<p) e~ 2iNv (hp 


( 68 ) 


is the matrix element of the disorder potential U(tp) 
mixing the states TV and —TV. In the effective model 
of Eq. (1671) we still neglect the matrix element corre¬ 
sponding to optical transition between the levels —TV and 
—TV — 1 by keeping in mind that such a transition is non¬ 
resonant for circularly polarized light. 

It is instructive to diagonalize the effective Hamilto¬ 
nian (1671) with respect to disorder potential. The corre¬ 
sponding basis states are, then, conveniently numerated 
by the index L or R and by n = TV or TV + 1 , 


T' 



VMi + \U 2 )’ 

E l (R) 

n 


" VMl + \Zn \ 2 )’ 

= D n ± ^Dl + \U n \ 2 , 


(69a) 

(69b) 


where the following notations are introduced 


D n 


-Era 


— <pne Oi 


U n 


D n + \f\U n \ 2 + D ^ 


(70a) 

(70b) 


For U n = 0, the functions 4'^ and correspond to 
right- and left-moving states, exp (irup) and exp (—imp), 
respectively. 

Using the basis functions e _ “ t 4'^ +1 , 4/^, and 
e -iuit x pi we re write the Hamiltonian given by Eq. (1671) 
in manifestly time-independent form, 


( E~ N V 0 V? N+1 \ 

V E^-U VC N 0 

o vi N e+ vi N i * N+1 

\U£jv+i 0 Vi^iN+i Ejf +1 -ui / 


( 71 ) 
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FIG. 7: Resonant (w w Ajv) optical transitions between levels 
of right- and left- moving electrons. Transitions LL, RL, and 
LR are induced by disorder. 


where we introduce 


fiab 


2 Vab sign(£ ab ) 

|^ab| + V ^'ab + 4| | 2 


S ah = W - A 


ab• 


(76) 


The result of Eq. (169aD can now be substituted into 
Eqs. ©, (|75aD . and (I75bl) in order to calculate the 
radiation-induced current. The calculation in the adia¬ 
batic case is very similar to those presented in Sec. IIII Al 
For disorder-induced splitting that is small compared 
to temperature one finds that the resonant radiation- 
induced current for w « Ajv is still given by Eq. (l36l) 
with 


Jn = Io ^ 

a,b 


I Pab 


i + |/W 


[(TV + 1/2 )A ab + B ab ! 2], (77) 


where we introduced 


where 


V = 


cue o 


2 ^ 


|6v| 2 )(l + |6v+i| 2 ) 


(72) 


In the absence of disorder one has 6v = 6v+i = 0; 
hence the only optical transition allowed is the RR tran¬ 
sition between the states TV + 1 and TV, which both cor¬ 
respond to right-moving electrons. The corresponding 
matrix element equals aeo/2 = e£(]R/2. 

In the presence of weak short-range disorder the origi¬ 
nal states for right- and left-moving electrons are mixed. 
As a result the RR transition occurs between the states 
^n+i and with the corresponding matrix element 
Vrr — V. In addition three more transitions emerge 
an d that are labeled 

as LR, RL, and LL, respectively (see Fig. [7]). As can 
be seen from Eq. ED the matrix elements corresponding 
to these transitions are Vlr = F6v+i, Vrl = Vt;* N , and 
Vll = V£h£n+ i- The corresponding resonant frequen¬ 
cies are given by 


A n h — E/y+i — E 


J AT> 


(73) 


where 2£“ is given by Eq. (I69bl) and the indices a, b take 
on R, L. For small V all four optical transitions are well 
resolved. If radiation frequency is close to one of the 
resonant frequencies A ab one may again use a two-level 
approximation described by the effective Hamiltonian 


K* b 


f zpa 

H ab = ( N b 

. Vab L n+1 


— UJ 


(74) 


which acts in the space spanned by the functions ’I' 1 ^ and 
t '&N+i- The eigenfunctions of the Hamiltonian (1711) 


are readily found as 


vl/ a - 
'a n ~ 


^jv - Pabe~ iut * b N+1 


\J 1 + |/?a 




N+l 


Pab^N+E ^AT+1 


\/l + \Pa 


Arr = —An = Brl = —Blr = Xn+i — Xn, 
Brr = —Bll = Arl = —Alr = Xn + Ajv+i, 


and 


An — 


1-16, 

i + 16 , 


12 * 


(78a) 

(78b) 


(79) 


The derivation of Eq. ED. which describes 4 resonances 
at frequencies A a b, has been based on several important 
assumptions. 

First of all the 4-level approximation used to justify 
Eq. (1671) is valid if disorder mixes nearly degenerate levels 
n and —n only, i.e.,those which have opposite chiralities 
in the absence of disorder potential. These levels are sep¬ 
arated by energy D n defined in Eq. (1 70 all , hence the mix¬ 
ing is controlled by the parameters U n /D n = Un/rupEo 
for n = TV, TV + 1. The admixture of other levels is weak 
as far as U n -C ueq, which is the central condition for 
the validity of Eq. ED- Note, however, that the relation 
between U n and D n can be arbitrary. 

It has been also implied that the resonance frequencies 
A ab arising due to the splitting of the TVth resonance of 
the clean ring do not overlap with the frequencies arising 
from the splitting of the (TV+ l)th resonance. This yields 
the condition Ncf) -C 1. 

Finally, we assumed that the radiation is sufficiently 
weak so that all four resonances predicted by Eq. ED 
are well separated. For weak disorder U n -C D. n the 
latter requrement is satisfied if a -C <j). 

The structure of resonances, which follows from 
Eq. (1771) . is shown in Fig. © assuming the limit of weak 
disorder |£„| <C 1 (or X n ~ 1). In this limit one can ne¬ 
glect the disorder-induced level repulsion hence the reso¬ 
nance frequencies are set by 


(75a) 

Arr = 

(75b) 

Alr = 
Arl = 




A ll = A N 
n + 2(TV + l)eo<6 


■ 2e n 


(80a) 

(80b) 

(80c) 
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The width and height of each of the resonances are deter¬ 
mined by the corresponding matrix element V a b. The cur¬ 
rent direction at resonance is also different. Remarkably, 
the amplitudes of LR and RL resonances are enhanced 
by a factor N « iip as compared to those of RR and LL 
resonances provided weak disorder regime |£ n | -C l/N. 
In this regime the LR and LL transitions correspond to 
antiresonances, i.e,the direction of current is opposite to 
those at RL and RR resonances. Thus, weak short-range 
disorder does not suppress or smoothen the RR resonant 
peak in current but leads instead to the appearance of 
two sharp and intense resonances of opposite chirality 
corresponding to LR and RL transitions. 

For stronger but still sufficiently weak disorder, such 
that l/N -C |£ n | -C 1, the amplitudes and signs of the LR 
and RL resonant peaks do not change, while the RR and 
LL resonances are strongly enhanced. The correspond¬ 
ing amplitudes are of opposite sign and proportional to 
i-^GCivI 2 — |£iv+i| 2 )- The absolute sign of the factor 
|£jv | 2 — |£/v+i| 2 depends, however, on disorder realiza¬ 
tion and cannot be predicted. Thus, resonant optical 
excitations in the ring can be used to probe mesoscopic 
fluctuations of disorder. More specifically, the sign of 
RR resonance might change for different disorder real¬ 
izations. 

Finally, we notice that the amplitudes of all 4 resonant 
peaks decreases provided disorder becomes so strong that 
U n D n and, consequently, A,,, -C 1. 

With increasing radiation intensity individual reso¬ 
nances start to overlap. For a </> only RR and LL 
resonances overlap while LR and RL resonances remain 
well separated. For a N<p all 4 peaks overlap and form 
a wide resonance which in the first approximation is de¬ 
scribed by Eq. (1551) . The effect of weak disorder remains 
small since the levels n and —n are no longer degenerate 
even at cf) = 0 due to radiation-induced level repulsion. In 
this case, the effect of disorder can be taken into account 
by the standard perturbative analysis using Eqs. (132a!) 
and (132bl) as zero approximation. Calculating the per¬ 
turbative corrections up to second order with respect to 
disorder potential we obtain 


Jn = 


a 


\J' ol 2 + Sfj(\ <5/v | + o? + 
8A/(|£/jv| 2 - \Un+i\ 2 ) , 

e o(IM + \[® r 


S 2 n ) 2 


• ( 81 ) 


This result is valid for a max{N<f>, '/N\U n \/eo}. In¬ 
terestingly, the sign of disorder-induced correction to cur¬ 
rent also depends on particular realization of random po¬ 
tential. For the case \Un\ < |bGv+i| the resulting current 
is plotted schematically in the Fig. |as a function of 
frequency. 

The case of nearly half-integer flux piercing the ring 
may be considered in a similar fashion. Assuming that 
(f> ss 1/2 we find that the energy levels are also double de¬ 
generate, but the disorder potential U(<p) mixes the level 


Ir ad 



FIG. 8: Disorder-induced splitting of N —th resonance into 
four peaks. Amplitudes of RL and LR peaks are enhanced 
by a factor N ~ uf. 


/rad 



FIG. 9: Structure of A—th resonance for Ncj> -C a. 


N , which has a positive chirality, with the level — (7V+1), 
which has a negative chirality. Corresponding matrix el¬ 
ements differ slightly from those given by Eq. (l68l) , 

Un = 2-J U(<p)e~ 2i( > N+1 ^d<p, (82) 

which does not make, however, a difference for the struc¬ 
ture of resonances. Still, the width, height and position 
of resonant peaks change accordingly. 


B. Long-range disorder 

In this subsection we consider the effect of static long- 
range disorder U(ip) that does not lead to scattering 
between left- and right- moving electron states. Con¬ 
sequently it is still convenient to analyze the effect in 
the rotating reference frame. The corresponding elec¬ 
tron wave function y(0) yields the Schrodinger equation 
m which is equivalent to those for a quantum physical 
pendulum. In contrast to the laboratory frame, the cor¬ 
responding potential in the rotating frame is no longer 
static but oscillates with the frequency w, 

OO 

[%) = U(6 + Ut) = Yl Un e in( - e+ut) . (83) 

n =—00 

Since the potential is smooth, i.e.,it does not change es¬ 
sentially on the scales of the order of the electron wave¬ 
length, the quasiclassical analysis is justifiable. The 
problem is, therefore, reduced to that of a classical physi¬ 
cal pendulum subject to a fast-oscillating potential. Such 
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FIG. 10: Levels captured in the chaotic layer (shown in gray) 
are randomly populated and depopulated leading to current 
fluctuations. 


a model is often considered in textbooks as the simplest 
example of a system that shows chaotic behavior. 44 

It is well known that the effect of the fast-oscillating 
potential is negligibly small everywhere except for a nar¬ 
row strip in the phase space that ” dresses” the separatrix 
(a curve which separates oscillating and rotating pen¬ 
dulum states). Such a strip is called the chaotic layer. 
Within the chaotic layer the physical pendulum jumps 
randomly between dynamically localized and delocalized 
trajectories (in the phase space) thus showing a chaotic 
behavior. The width of the chaotic layer can be estimated 

44 

as 2 s 


T ch oc |f7i| exp > ( 84 ) 

where U\ is the amplitude of the first Fourier harmonic 
of the oscillating potential [see Eq. ©] and is the 
characteristic energy scale corresponding to pendulum 
frequency at the point of equilibrium. 

The analogy between a quantum disordered ring sub¬ 
jected to circularly polarized light and the physical pen¬ 
dulum in oscillating potential suggests that the long- 
range disorder may play a role only in the limit of strong 
coupling to electromagnetic field, i.e., for a 1. In this 
case the number of levels perturbed by the radiation- 
induced potential W is proportional to y/a 3> 1. This 
justifies the quasiclassical approach which suggests that 
f l = y/a£o is the characteristic energy scale which enters 
Eq. JH3. 

The chaotic layer separates the regions with localized 
and delocalized states. For strong coupling the chaotic 
layer is confined to an energy interval around e = 0 of the 
width 5e ~ T c h (see Fig. [2j. The energy levels captured 
by the chaotic region correspond to the values of n such 
that \n — n w | ~ yfa as shown in Fig. [10] The number of 
levels within the chaotic layer can be estimated as 


(85) 


where p c h is a density of electron states within the 
chaotic layer. Simple quasiclassical analysis of Eqs. m 
and (fl5l) shows that p(e) ~ ln(a/e)/eo yfa. in a vicin¬ 
ity of the separatrix, so that p c h can be estimated as 
Pch ^ ln(a£o/r ch )/e 0 y/a. 

The chaotic behavior leads to fluctuations of dc cur¬ 
rent due to random jumps within the chaotic layer. The 
amplitude of such fluctuations is given by 

A I ~ (86) 

an 

In order to estimate the derivative dl n /dn we take ad¬ 
vantage of the perturbative result given by Eq. (l27l) which 
is taken at the boundary of its applicability range, i.e. for 
\n — n^l = yfa. In this way we find dl n /dn ~ Iq. As¬ 
suming that r ch <C aeo we obtain from Eqs. (1841186ft that 


* t r r ch ln(ae 0 /r c h) 

~ l 0 - 7 =- 

£0 V a 


(87) 


The fluctuations of current due to chaotic dynamics re¬ 
veal themselves as a burst noise, i.e. a random telegraph 
signal. An electron captured within the chaotic layer 
makes a random hop during a time p c h which is a char¬ 
acteristic time required for an electron to travel along a 
typical trajectory within the chaotic layer. Sudden jumps 
in the current of a magnitude A I correspond to a hop¬ 
ping event for any of 5n available electrons. Thus, the 
characteristic rate of current jumps can be estimated as 
l/r c h « Sn/p c h- Hence r c h — l/r c h as can be seen from 
Eq. (1551) . 


VI. BALLISTIC MULTI-CHANNEL RING 

The generalization of our results to the case of ballis¬ 
tic multi-channel rings is rather straightforward. In the 
absence of radiation the energy levels in a multi-channel 
ring are given by 

eL°1 = £ oy+<r, (88) 

where the additional index m numerates the subbands 
(channels) due to transverse quantization. The specific 
expression for subband energies e™ depends on the con¬ 
fining electrostatic potential in the transverse direction. 

The contribution of the mth subband to optically in¬ 
duced dc current can be easily found from equations 
obtained in the previous sections by replacing the Fermi 


energy with the energy 

E% = E f - e?, (89) 

and by redefining the parameters such that 

uf -a rip = E™/£q, (90a) 

Aj, -a A£ = e 0 <, (90b) 

Sui -A 6oj m = T/n%. (90c) 


6n — PchTch} 
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For weak coupling a -C 1 the position of resonances 
would correspond to the energies E^ F1 m -E<®. For 
each channel m the transverse energy drops out from this 
expression hence the parameters S n entering Eq. (E3 are 
still given by Eq. (129al) . The existence of several conduct¬ 
ing channels would affect, however, the envelope function 
which is given by Eq. (1391) for a single-channel ring. This 
result is, however, readily generalized with the help of 
the substitutions Performing the summation over 

all subbands we obtain the current 


^env(^) 


t rnF A m 

1 0 _ 

8T “'o - A£?)/2(5u; m ] ’ 


( 91 ) 


where mj? is found from the equation = Ep. 

The resonant peaks in the envelope function I env (uj), 
which correspond to the mth and (m + l)th subbands, 
are well separated as far as A ™ +1 — A™ Su> m . The 
latter inequality is equivalent to the condition 


r)p m 

tt » r - < 92 > 

which can be obtained with the help of Eqs. (liiOl) using 
that A™ +1 — A™ Ri <9A™/<9 to. Thus, the subbands pro¬ 
vide well separated resonant contributions to the current 
provided the temperature is small compared to the inter¬ 
subband spacing. 

For higher temperatures inequality (1921) is violated and 
the resonant peaks in the envelop function I env (uj) start 
to overlap. Finally, at very large temperatures, such that 
T de™/dm for any m < mp, all resonances merge in 
a single wide peak which can be described by Eq. m, 
where the summation over m is replaced with integration, 


fenv(^) 


- / 

8T J 


dm 


A m 
^F 


cosh 2 [(w — A£l)/2(5u; m ] 


(93) 


The integral is readily estimated for the simplest model 
of the ring of a finite width a assuming that the effect 
of the confining potential is properly accounted by using 
periodic boundary conditions in the transversal direction. 
For such a model we obtain 


e ±= e ^ (94) 

where ex = (2irh) 2 / M a 2 . We also find 


E™ 

= Ep(l — m 2 /m 2 F ), 

(95a) 

A^ 

= A F ^l~m 2 /m 2 F , 

(95b) 

5uj m 

= 8oj/\Jl — m 2 /m F , 

(95c) 


mF = \j2Ep/£p. 

(96) 


Using that Ap Suj (this inequality is fulfilled for Ep 
T) one can find the asymptotic behavior of the integral 
in Eq. (l93ll as 


fenv (^ ) 


J < 


Ioy/T/e± I 

4 n* j 2uiy/s^ /A~ B 


V \Z A f~ u 


uj — A p 3> Suj, 

|w — Ax| <C Suj, 

A]?-w> Suj, 

(97) 


where C = / 0 °° dx/ cosh 2 (a; 2 ) ft 0.95. 

The result of Eq. (1U71) predicts exponential decay of 
current for large frequencies, w - Aj? > Suj. One can, 
however, see that such behavior is limited by uj < A p + 
5uiIii(Ef/T). For larger values of uj, I env decays in a 
slower power-law way. 

Comparing Eq. (1971) with Eq. (l39l) we find that at high 
temperatures the optical response in the multichannel 
ring is enhanced by a factor ^JT /ex- 


VII. DISCUSSION AND CONCLUSION 


In this section we estimate the value of the current 
I ra d induced in a semiconducting nanoring by circularly 
polarized light. We also discuss related problems which 
are to be addressed in future, and summarize the results 
obtained. 

The only parameter of the theory which depends on 
specific material properties is the electron mass M. To 
make estimates we take the standard value of the effec¬ 
tive mass for GaAs, M ft 0.07mo, where mo is the free 
electron mass. We also take the radius of the ring to 
be R = 100 nm. For a single-channel ring we, then, 
obtain using these parameters that £o ~ 10 -4 eV and 
Iq ft 3.7 x 10 _ 9 A. For a wide range of Fermi ener¬ 
gies Ep = 0.02 — 2eV, we find uf ft 20 — 200 and 
Ax ft (0.2 — 2) x 10 -2 eV. We see that inequality Eq. J3]) 
can be easily satisfied for not too large temperatures. The 
corresponding resonance frequency turns out to be in the 
THz range, / = w/27t = Ap/2ir ft 0.5 — 5 THz. The cou¬ 
pling to electromagnetic field becomes stronger for suffi¬ 
ciently low fields such that a = 1 for £q ft 10 V/cm. The 
maximal value of the current in the weak coupling regime 
(this value is reached in LR and RL disorder-induced 
resonances) is estimated as Iotif ft (0.8 — 8 ) x 10“' A. 
In the very strong coupling regime, a > n|, the cur¬ 
rent increases to reach a maximal value / max ft Ion 2 F ft 
(1.5 — 150) x 10 ~ 6 A. The conditions for the very strong 
coupling regime are satisfied only for sufficiently strong 
fields: £ 0 > (4 — 400) x 10 3 V/cm. 

Let us now estimate the magnetic field induced by the 
current / ra d- For a single ring, the field in the center of 
the ring is given by a simple formula 


B = 


27r/ ra d 
cR ’ 


(98) 


where c is the speed of light in the vacuum. In the weak- 
coupling regime, / ra d ft IoUf, we estimate B ft (0.5—5) x 


where 
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10 _6 T. In the strong-coupling regime the maximal field 
B max ~ (0.1 — 1) x 10- 3 T is reached for J rad ~ / max . This 
field can increase further in a multi-channel ring and/or 
by using three-dimensional arrays of rings. Another way 
to increase the effective magnetic field generated by the 
ring is to make the Fermi energy and, consequently, the 
parameter uf larger. 

In our analysis we focused on the dependence of the 
J ra d on the frequency u of incoming radiation and found 
that circular current might show sharp resonances. Im¬ 
portantly, the dependence of / rad on the magnetic flux 
(f> also reveals sharp peaks for a given ui provided inter¬ 
action constant a is sufficiently small. Indeed, the posi¬ 
tions of resonances, which are shown in the upper panel 
of Fig. [5] depend on magnetic flux. Increasing magnetic 
flux by the flux quantum, (/—>■</> + 1, is equivalent to 
the substitution S n —> 5 n -1 in Eq. ED that describes 
the dc photoresponse for the case of adiabatic radiation 
switching. (The same remains true for the ring coupled 
to a thermal bath provided weak coupling to electromag¬ 
netic radiation.) Thus, the current / rad is a periodic 
function of <j> with the period 1 as expected. It is ev¬ 
ident from the consideration above that there exists a 
single sharp peak in </> dependence of / rad in the interval 
0 < (j) < 1. The ratio of the maximal value of the cur¬ 
rent, /™ d x in this interval to the flux-averaged current 
(Aad )<t> is as large as 1 /a and 1/y/a 2 + F 2 /2 for the case 
of adiabatic radiation switching and thermal bath cou¬ 
pling, respectively. This ratio also gives an estimate for 
the number of harmonics that are effectively contributing 
to the Fourier expansion of the radiation-induced cur¬ 
rent / ra d = exp(27 Tim<j)). The bigger the ratio the 

larger the number of relevant harmonics with a large am¬ 
plitude that can be observed in experiment. 

The dependence of the circular current on (j) is essen¬ 
tially different in the strong-coupling regime, a> 1. In 
this case, all harmonics I m for m ^ 1 are small com¬ 
pared to Iq. In particular, I m /Io oc exp(—27r 2 |m|u*) or 
Im/Io oc exp(—27r|m|F) for adiabatic radiation switch¬ 
ing (assuming 1 <C n* -C y/a) and for the case of ther¬ 
mal bath coupling (assuming T = T^,, and 1 < a < 
T 2 ), respectively^ Hence, in the case of strong cou¬ 
pling to the radiation the response is given by a large 
flux-independent quasiclassical contribution and a small 
quantum correction oscillating with <f>. The latter is dom¬ 
inated by the contribution of harmonics with m = ±1. 

In the presented analysis we ignored the effects of the 
electron-electron interactions that may not be negligible. 
Even thought the detailed study of interaction-induced 
effects is a complex task that falls outside the scope of 
the current paper some qualitative predictions can be 
already made. 

One may expect that in a single-channel ring at suf¬ 
ficiently low temperatures the main effect of electron- 
electron interactions is to renormalize the value of the 
coupling strength a. Such renormalization will likely re¬ 
sult in the suppression of a by a factor ( T/Ep) 9 , which 
is characteristic for the Luttinger liquid behavior, where 


g is a dimensionless interaction constant. This effect can 
be taken into account by replacing the coupling constant 
a with its renormalized value that would not change es¬ 
sentially the predictions of our analysis. 

If interactions are sufficiently strong, less trivial effects, 
which are related to the charge quantization in a finite 
geometry, may show up. As was first demonstrated in 
Ref.the electron-electron interactions in a ID ring give 
rise to an effective contribution to magnetic flux that is 
proportional to both the interaction constant g and the 
imbalance Np—Np, where Np(Np) is the total number of 
right-(left-) moving electrons in an ideal ring. As a result 
sufficiently strong interactions in a clean system would 
lead to further splitting of the four resonances described 
in Sec. |V] 

A completely different but sizable effect of interactions 
is expected in multi channel rings for sufficiently high 
temperatures such that the electron-electron collisions 
dominate. This case is generally referred to as the hydro- 
dynamic regime. In this regime elementary excitations in 
the ring are dominated by plasmons. The corresponding 
plasmonic resonance in the dc current has a width which 
is much smaller than the resonance width in the ballistic 
noninteracting ring studied above. The decrease of the 
linewidth is due to the motional line narrowing caused 
by intense electron-electron collisions. Such and other 
interaction-related phenomena will be studied elsewhere. 

To conclude, we developed a theory of the inverse reso¬ 
nant Faraday effect in quantum rings. We demonstrated 
that a circularly polarized radiation with the frequency 
ui induces a dissipationless dc current / rad in a quantum 
ring pierced by magnetic flux <f>. The current yields the 
symmetry, Eq. ©, so that the direction of the optically- 
induced current is sensitive to helieity of the incoming 
radiation. 

For the case of weak coupling the current J ra d(w, (f>) re¬ 
veals sharp resonances as a function of uj for a given flux 
</>. These resonances can also be observed by changing 
the flux for a fixed frequency of light w. Analytical ex¬ 
pressions for the radiation-induced current are obtained 
for two different cases: (i) an isolated ring under the as¬ 
sumption of adiabatic switching of light intensity, and (ii) 
a quantum ring weakly coupled to the thermal bath. 

The nonresonant current is found to be proportional 
to the squared amplitude of light / rad oc £$ in agree¬ 
ment with the conventional theory of the nonresonant 
inverse Faraday effect. The current is, however, strongly 
enhanced at a resonance so that its maximal value does 
not depend on the intensity of light (in the regime when 
dissipation is negligible). 

For the case of strong coupling multiple resonances in 
/ ra d(w,</) merge into a wide peak with a width deter¬ 
mined by the spectral curvature. The amplitude of the 
peak increases with £q as £ 0 ' for small light intensity. It 
is proportional to £q for moderate intensities and finally 
saturates in the limit of hight intensity of light. The sat¬ 
urated value of the current scales as n 2 F with the total 
number of electrons in the ring uf- 
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Weak disorder is shown to affect the dependence of J ra d 
on frequency in a highly nontrivial way. In contrast to 
naive expectations a weak short-range disorder does not 
suppress resonances but leads instead to the appearance 
of additional resonant peaks of different polarity. These 
sharp resonant features are suppressed only by relatively 
strong disorder potential. Thus, we find that the inverse 
Faraday effect is generally very sensitive to the quality of 
the ring. 

The long-range disorder does not affect the picture for 
the case of weak coupling to light while it becomes essen¬ 
tial for the case of strong coupling. The main effect of 
long-range disorder is to induce a chaotic behavior of the 
system in the vicinity of the separatrix that divides the 
phase space into the regions with dynamically localized 
and delocalized states. The radiation-induced current 
Trad is shown to fluctuate due to random electron hop¬ 
ping within the narrow chaotic layer in the phase space 
of the system, which ” dresses” the separatrix. Such fluc¬ 
tuations lead to the burst noise in the optical dc response 
and power dissipation. 

Finally, we generalize some of the results obtained for 


the case of a multichannel ring. We demonstrate that 
at low temperatures the response originating in different 
propagation channels is well separated in frequency so 
that the spectrum of different subbands can be resolved 
in experiment. At higher temperatures the resonances 
overlap but the overall response is enhanced by a factor 
\jTje _l as compared to the case of a single-channel ring. 
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